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Abstract 

This paper contains a detailed derivation of the photoinduced current density at 
third order in the couphng between a semiconductor and a multifrequency photon 
field, starting from its standard textbook expression as a third order time integral of 
a triple commutator. Due to a major intrinsic problem linked to this triple commu- 
tator, such a derivation has been made possible quite recently only, thanks to the 
tools developed in the composite-boson many-body theory we have just constructed. 
The photoinduced current density is shown to ultimately read in a compact form, 
in terms of the Pauli and Coulomb scatterings for exciton-exciton interactions in- 
troduced in this theory. Representation in Shiva diagrams is also given to better 
grasp the physics of the various contributions. 

PACS number: 71.35.-y 
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1 Introduction 



Nonlinear effects induced by unabsorbed photons are Icnown to be of high technologi- 
cal interest due to the timescale on which ultrashort laser pulses allow them to operate. 
A quite usual way to describe these optical effects is through non-linear susceptibilities 
[1-17]. Since photons do not interact directly, these nonlinearities can only come from 
interactions between the matter excitations to which these unabsorbed photons are cou- 
pled. In the case of semiconductors, these matter excitations can be conveniently seen as 
virtual excitons. Interactions between excitons have two quite different origins: Coulomb 
interaction between the carriers of these excitons, but also and mainly Pauli exclusion be- 
tween these carriers. This Pauli exclusion remained for decades quite difficult to handle 
properly. This is why excitons have been commonly treated as elementary bosons with 
effective scatterings dressed by a certain amount of carrier exchanges [18,19]. 

Over the last few years, we have developed a many-body theory for composite bosons 
made of two fermions, like the excitons, which, as a main goal, has the exact treatment 
of fermion exchanges between composite bosons [20,21]. This theory, in a natural way, 
generates "Pauli scatterings" which describe carrier exchanges between two excitons in the 
absence of carrier interactions. By combining these 2x2 scatterings, it is actually possible 
to describe all carrier exchanges which can exist between excitons, these multiple 
exchanges being nicely visualized through the so-called "Shiva diagrams" [21,22]. Since 
these Pauli scatterings are dimensionless by construction, they in fact control all optical 
nonlinearities when the detuning increases, due to a bare dimensional argument: Coulomb 
scatterings associated to Coulomb interaction between carriers are energylike quantities, 
so that when they appear, they must have energylike denominators which can only be 
photon detunings. This makes Coulomb processes negligible in front of pure fermion 
exchanges when the detuning increases. And indeed, in essentially all nonlinear effects 
we have up to now studied [23-29], Coulomb interaction between excitons plays a minor 
role. We will see that this is also true in the present work, the dominant term in the 
final expression of the photoinduced current density, Eqs.(5.1,2), again being the Pauli 
scattering term. 

Textbooks usually give the third order response to a photon field as a third order 
integral over time of a triple commutator [1,8,13]. This algebraic expression, although 
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nicely compact, immediately raises a major problem since, when developed, this triple 
commutator makes appear the Hamiltonian acting on two-electron-hole-pair states. As the 
exact eigenstates of two electrons and two holes are not known except for highly simplified 
Hamiltonians [5,11], the exact calculation of these terms is not possible. It is however 
necessary to somehow control them since they contain volume linear contributions which 
have to be extracted in order to show that they cancel exactly, for nonlinear susceptibilities 
are intensive quantities. This " cancellation problem", which remained open for decades 
[2-4], has been successfully tackled recently [25,26] using the tools we have developed for 
interactions between composite excitons. In this framework, the calculation of the third 
order susceptibility turns out to be rather simple because, in it, only enter two interacting 
excitons, not as in many other problems we can now address. The formal elimination 
of these volume linear terms relies on the possibility to, in an exact way, pass e*^* over an 
exciton creation operator. This "trick", which is given in Eq.(3.5), turns out to be useful 
not only for calculating nonlinear susceptibilities but also for any problem involving time 
evolution of exciton states. 

Once these volume linear terms are eleminated, we are left with an expression of 
the nonlinear response to a photon field which reads in terms of Pauli scatterings for 
carrier exchanges between two excitons and Coulomb scatterings for carrier interactions, 
as physically reasonable since nonlinear optical effects in semiconductors are expected 
to come from interactions between the virtual excitons to which photons are coupled. 
This result was reported in ref. [26]. A somewhat more detailed version can be found in 
ref. [25]. However, the procedure we first used to obtain it was rather heavy: We calculated 
explicitly the eight terms of the triple commutator and carefully combined all the volume 
linear terms to show that they exactly cancel. In writing the extended version of this 
work, we have realized that, by calculating the first of these commutators explicitly, we 
are trivially left with Pauli and Coulomb terms only, so that, by doing so, we immediately 
get rid of all the nasty volume linear terms. It is this far simpler derivation that we here 
report in details. In view of its formal simplicity, it becomes easy to consider a photon field 
with multifrequency, as necessary to possibly describe four-wave-mixing experiments in 
which different laser beams would be used. In this work, we however restrict to photons all 
having the same circular polarization. This allows us to drop the spin indices and makes 
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all notations far lighter. Extension to photons with different elliptical polarizations, in 
order to cover all possible configurations, will be given in a forthcoming paper, the goal of 
the present work being to concentrate on the fundamental aspects of the new procedure 
we propose. 

Due to a bare dimensional argument, the dimensionless Pauli scatterings cannot ap- 
pear in boson-exciton effective Hamiltonians [18,19] whatever the bosonization procedure 
to generate them is. As these Pauli scatterings control all optical nonlinear effects at 
large detunings, it is clear that it is not possible, through these effective Hamiltonians, 
to properly describe optical nonlinearities in semiconductors, nor to derive the correct 
expression of nonlinear susceptibilities. 

Various groups [7,9,10,12], through rather different procedures, have managed to tackle 
the third order susceptibility while keeping the exciton composite nature exactly. Since in 
the third order susceptibility, only enters the interaction of two excitons, it is not necessary 
to have at hand a quite general A^-body formalism to possibly solve this two-body problem 
exactly. In these works, the third order response to a photon field is usually approached 
not through the standard textbook expression for nonlinear susceptibilities, but through 
its time derivative. This turns out to be quite wise as the so-called cancellation problem, 
which immediately appears when writing this susceptibility through a triple commutator, 
is then totally avoided. Comparison between the present approach and these various 
procedures, within the same set of notations in order to make an easy link between all of 
them, is clearly of interest as the obtained results can appear as rather different at first, 
even to a careful reader. However, this comparison, which implies to go into each of these 
different procedures rather in details, is definitely too long to be included here. This is 
why we will present it in an independent paper. 

The present paper is organized as follows. 

In section 2, we settle the problem. 

In section 3, we outline the difficulty the "brute force" calculation of the triple com- 
mutator raises and give the trick which allowed us to overcome it. 

In section 4, we show how we can, within a few lines, get rid of all volume linear terms 
and immediately reach the physically relevant Pauli and Coulomb terms. We then show 
how to get a compact form for these two terms. 
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In section 5, we discuss the final result and give its representation in Shiva diagrams. 
In section 6, we conclude. 

This paper also contains two appendices. In the first one, we come back to the ap- 
propriate way to describe optical nonlinearities in solid state physics, to once more stress 
that, although seen in many semiconductor textbooks and publications, it is not the pho- 
toinduced dipole density but the photoinduced current density which has to be used in 
this description. The second appendix contains a compact rederivation from scratch of 
the formal expression of operator mean value at third order in coupling, for completeness. 

2 The problem 

The standard way to approach optical nonlinearities is to calculate the time dependence of 
physically relevant operator mean values due to the wave function change induced by the 
coupling of the system at hand to a photon field, this field being introduced adiabatically 
from t = — oo. As explained in appendix A, the relevant operator in solid state physics 
is not the dipole density P(r) as used in atomic physics, but the current density J(r), 
because wave functions for solid state systems are not spatially localized as in the case of 
atoms but extended over the whole sample, due to crystal periodicity. 

Since the relevant excitation operators for problems dealing with semiconductors are 
the exciton creation operators i?^, i.e., operators such that 

H,,Bl\v) = E.n,Bl\v) , (2.1) 

where Hgc is the semiconductor Hamiltonian and \v) the electron-hole vacuum, the in- 
terband contribution to the current density, which is a one-electron operator, can be 
conveniently written in second quantization, in terms of these exciton operators, 

J(r) =Ej-W^m + h.c. . (2.2) 

m 

As rederived in Appendix A, the prefactor jm,(r) is related to the Kane vector G [30] 
for valence- conduct ion transitions, given by Eq.(A.8), through jm(i') = — G e***^™'''(r = 
Oli'm) L~^^^, where is the center-of-mass momentum of exciton m, v^n its relative 
motion index, L the sample size and D the space dimension. 
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The interaction Wt between a semiconductor and a photon field can also be writ- 
ten in terms of exciton operators. This actually is far more convenient than using free 
electron-hole-pair operators as commonly done, since excitons are the physically relevant 
excitations in semiconductors. For a multifrequency photon field with vector potential 
A(r,t) introduced adiabatically from t = — oo over a timescale 1/e, 

A(r, t) = e'' e'^'^^'-^^-"^ + c.c. , (2.3) 

j 

where {uj, Qj) are the frequency and momentum of the photons at hand, this interaction 
reduces, if we only keep resonant terms (which corresponds to the so-called "rotating wave 
approximation" [31]), to 

Wt = Ut + Ul , (2.4) 
where Ut destroys an exciton according to 

Ut = Y.^n{t)Bn, (2.5) 

n 

with a time- dependent weight given by 

x„(t) = e^*5:/i(^V-^*. (2.6) 

i 

The prefactor ^^^^ in Xnif) is related to the same Kane vector for valence-conduction 
transitions appearing in the current density, through 

/x(f) = G.A,.5Q„,Q^(r = 0|z/„,)L^/^ (2.7) 

Due to the (r = Oli^^) factor, photons are only coupled to excitons with S symmetry, 
for their wave function to differ from zero at r = 0, the largest coupling being with the 
ground state because its wave function at the origin is the largest. 

The interaction expansion of the current density mean value (J(r))t = J(r)|?/'j) for 
a system wave function \'il)t) which obeys the Schrodinger equation i§i\i^t) = {Hsc+Wt)\ilJt) 
only has odd order terms in the semiconductor-photon interaction, since both J(r) and 
Wt create or destroy one exciton. For initial state \il)t^_^) taken as the electron-hole pair 
vacuum the third order term in Wt is known [1,8,13] to be given by the integral of 
a triple commutator (see Appendix B). Due to Eqs.(2.1) and (B.5), (J*'^-*(r))f ends by 
reading as 

(J^'Hr))* = Eb"^W7m(t)+j;(r)7;(t)] , (2.8) 

m 
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where the prefactor 7j^(t) is given by 

l*mit) = H? f dh r dt2 r dhSmih,t2,h;t) . (2.9) 

J—oo J —oo J —oo 

Sm{ti-, h^t^] t) has a nicely compact form, 

SUti,t2,h;t) = {v\ [[[Blit),W{t^)],W{t2)],Wih)] \v) , (2.10) 

which, however, is not nice at all for analytical calculation, as shown below. In the above 
equation, Bj^(t) and W(t) are the Heisenberg representations of operators Bj^ and Wt, 
namely, Z{t) = e^^sct^^-^Hsct _ (^^Q^e that the prefactors of j*^(r) and jm(i") in Eq.(2.8) 
must be complex conjugate since (J'^^)(r))t is a real quantity). 

The problem is to get a compact expression for 7m (^)- In the next section, we are 
going to show why the "brute force" calculation of the triple commutator appearing in 
Sm{tii t2, t^, t) raises a major technical problem which let it open for decades. Thanks to 
the many-body theory for composite-bosons we constructed [20,21], we were able to over- 
come this difficulty [25,26] and to derive a compact expression of {3^^\v))t in terms of the 
Pauli and Coulomb scatterings of this theory. While the expression of (J*^^^(r))t we then 
obtained is fully correct, it turns out that, in preparing the present manuscript, we have 
found a way to greatly simplify our first calculation. This more elegant derivation is given 
in section 4. It however is of interest to put this new calculation in its "historical" context 
by briefiy presenting the "brute force" calculation of the triple commutator, in order to 
grasp the fundamental difficulty raised by a naive approach to third order susceptibility. 



3 Brute force calculation 

The brute force calculation of 7m (^) corresponds to expand the triple commutator of 
Eq.(2.10) and to calculate its eight terms separately. This is actually what we first did 
and reported in refs. [25,26]. This expansion makes appear "easy terms" and "tricky 
terms" . Easy terms are trivial to calculate because their intermediate state is the vacuum. 
Equation (2.1) leads to 

(H^„,(r4)St3(r3)5„,(r2)5t^(ri)|i;) = {v\BUT,)Bl{n)\v) {v\B^,{T2)Bl{r{)\v) 
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In contrast, in tricky terms, the intermediate state is a two-electron-hole-pair state. 



X 



(i;|i?„,i?„3e-^-(^^--)«K;) , (3.2) 



so that these tricky terms are impossible to calculate exactly as the semiconductor eigen- 
states for two pairs are unknown. 

The trouble is that the easy terms generate volume linear contributions which cannot 
exist in the final expression of the linear susceptibility, as this quantity is intensive. The 
volume linear contributions of the easy terms thus have to cancel out exactly with similar 
ones coming from the tricky terms. In order to prove it and extract these volume linear 
terms from the tricky terms, it is necessary to find a way to manipulate these tricky terms. 
Up to our many-body theory for composite bosons and the tools it provides, this was not 
possible; this is why the calculation of this triple commutator remained an open problem 
[2-4], except for very simple Hamiltonians for which the two-pair eigenstates could be 
analytically determined [5,11]. Let us now outline how the volume linear part of the 
tricky terms can be extracted. 

In the RHS of Eq.(3.2), we can commute e"*^""^ and -B^^ by using the integral repre- 
sentation of the exponential, namely. 



hoo 



_ J -i(x+ir]tO+)t 

J-oo [-2mr]t) X + ir]tO+ - Hsc 
with rjt = sign(t), and the key equation for correlations with excitons [20,21], namely, 

^ Bi = (bI + VA ^ J ^ , (3.4) 



O — Hsc V a — Hsc J a — Hsc — En 

being the creation potential for Coulomb interaction with the exciton n, defined as 
[HscBl] = EnBl + V^. This leads to (cf. ref. [21] section 6.2) 

^-^Hs.t^]^ = Bl e-^(^-+^")* + V^{t) , (3.5) 

^t(i) _ /"^^ e-^(^+^''*°+)* 1 ■ (3 6) 

" J-oo {-2mr]t) X + ir]tO+ - Hsc x + ir]S+ - Hsc- E^' 

When used in Eq.(3.2), Eq.(3.5) generates two terms. The first one, which now has Hsc 
acting on one exciton, readily gives, due to Eq.(2.1), 



The above scalar product can then be expressed in terms of Pauh scatterings for carrier 
exchanges between excitons as [20,21] 

{v\Br.,Br..^Bi^Bl\v) = 6n,,n, " A n?)) + (^4 ^3) . (3.8) 

If we then collect all the 6 parts of the various tricky terms produced by this procedure, 
it is possible to show that they cancel out all the easy terms. We thus are left with 
contributions from the tricky terms coming from the A part of Eq.(3.8) and from the 
V^{t) term of Eq.(3.5). This makes the physical origin of these remaining terms quite 
clear: 7m (^) has one part coming from Pauli exclusion and one part coming from Coulomb 
interaction. 

It is possible to calculate the Pauli contribution to 7m (^) analytically in terms of 
exciton energies and Pauli scatterings A (^"J with exciton m. In contrast, the fact that 
two-pair states play a role in this problem somehow remains, since the Coulomb part of 
7m(t) contains scalar products like (f |i?„^i?„2 ^ _^ V}^Bl^J\v). By using one of the 
key equations of the composite-boson many-body theory [20,21], namely, 

[VLBi]=Y.i(^,l)BlBl, (3.9) 

where ^ ^ is the direct Coulomb scattering of this theory, we can rewrite this scalar 
product as 

{v\Bn.B^, ——^ — ViBllv) = E e C; 2) {v\Bn,B^ ^ Bl Bljv) . 

(3.10) 

The exact calculation of the above matrix element is not possible because it imposes 
the knowledge of the whole two-pair eigenstate spectrum analytically. We can however 
estimate such a term in two limits, 

(i) close to a biexciton resonance: If |XX) is the biexciton ground state with energy Exx, 
Eq.(3.10) then reduces to 

{v\B^.B^,—^^ VlBl\v)= E , , "'^"^^^ {v\B^MXX){XX\Bi,BUv)-, 



(3.11) 

(ii) at large detuning: Equations (3.4), (3.8) and (3.9) allow to expand Eq.(3.10) in ratios 
of Coulomb scattering divided by detuning which are small in the large detuning limit. 
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Even if the above procedure is rather straightforward once Eq.(3.5) is given, to indeed 
show that all the volume linear terms cancel out exactly relies on an analytical calculation 
which turns out to be extremely heavy. It is then reasonable to think that if such a 
cancellation has to exist, there must be a clever way to show it, which avoids going 
through the calculation of these eight terms explicitly. 

This is what we present in the next section. Due to the formal simplicity of this 
calculation, it becomes easy to include a multifrequency photon field, the result we first 
reported being obtained for one photon frequency only. 

4 Better procedure 

We come back to Eq.(2.10) and first calculate the commutator [Bl^{t),W(ti)]. This is 
easy to do thanks again to Eq.(3.5). We find that this commutator contains a Coulomb 
term in which enters the creation potential V^, a Pauli term in which enters the deviation- 
from-boson operator Dnm. and a scalar. Since the scalar disappears when included in the 
triple commutator of Eq.(2.10), this readily shows that 7m(^), as well as all higher order 
susceptibilities which read in terms of higher order commutators, only have Pauli and 
Coulomb contributions. This result, obtained within a few lines, is already quite nice 
for the understanding of the physics which controls optical nonlinearities. Let us first 
calculate this commutator. 

4.1 Explicit calculation of the first commutator 

By using the definition of Heisenberg operators, this first commutator [Bj^{t),W(ti)] = 
Rm{ti',t) expands as 

To calculate it, we pass e~^^scit-ti) ^-y^^. ^^i^ ^^^^ term, and e-^^sdh-t) ^^^^ 

the second term. According to Eq.(3.5), this leads to split Rm(ti;t) as 

Rm{tu t) = RZ°"\ti; t) + e*^™(*-*i)e*^-*^ [Bl, Wt,]e-'"^^'' . (4.2) 
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The first term is a Coulomb term as it contains the creation potential V^. Its precise 
expression reads 

The commutator in the second term of Eq.(4.2) can be calculated by using one of the key 
equations for composite-boson many-body effects [20,21], namely, = Sn,m — Dnm- 

From the definition of Wt (Eqs. (2.4,5)), we find, since [Bj^, [//j = 0, 



ni 

= -Xmih) +'^Xn^{ti)Dn^rn ■ (4.4) 
ni 

By inserting this equation into Eq.(4.2), we end with 

Rmih; t) = R'rih; t) + Rl:-\h; t) + Rg"-\h; t) , (4.5) 
where R^^''(ti;t) = -Xm(ti)e*'^'"(*"*i^ is a bare scalar, while i?^'^"^'(ti; t), given by 

= e^^-*^ 5:x„,(ti)Ann^e-^^-*^e-*^'"(*^-*) , (4.6) 



ni 



is a Pauli contribution since it reads in terms of the deviation-from-boson operators Dn^rn- 
When inserted into Bl^if), W{ti) , W{t2) , the free part R^:^\ti; t) which is a scalar 
disappears; we are thus only left with Pauli and Coulomb contributions to 7m(i), as 
obtained (after a lot of efforts) through the brute force approach. Note that the same 
commutator Rm,{ti]t) actually shows that, not only the third order susceptibility, but all 
higher order nonlinear susceptibilities are only made of Pauli contribution and Coulomb 
contribution, i.e., they only come from interactions with exciton m through Pauli exclusion 
and Coulomb interaction, as derived from very fundamental physical arguments based on 
the fact that photons cannot directly interact. 

4.2 Formal structure of 7m(t) 

Due to Eq.(4.5), 7;;(t), splits as 7;;^''^''(^) +lm°''Kt), which are obtained from Eq.(2.9) 
with Sm{ti,t2,t3;t) respectively replaced by 



S^''''%,t2,h;t) = {v\ Rl'r\h;t),W{t2) ,W{h) \v) (4.7) 



)Pauli 
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and by S^°''\ti,t2,t3;t), which reads as ^^^"^'(ti, ts, ^s; with R^''^'{ti;t) replaced by 

Since [[Ai, A2] , A3] = {AiA2-A2Ai)A3-A3{AiA2-A2Ai), each of these ^^(^1, ^2, ^3; t) 
contains four terms. Fortunately, most of them reduce to zero. Indeed, from D„/„|f ) = 
which is the key property of the deviation-from-boson operator [21], along with -D^/„ = 
Dnn', we readily see that {v\R^^^\ti] t) and R^^^^ti, t)\v) reduce to zero. Since Ut\v) = 0, 
we are thus left with 

-S^^^\h,t2,t3;t) = (t;|f/,,e-^^-*^<r"(^i,^)e^''^^*^t//3l^) + (^2 - h) , (4.8) 

where R^'^^\ti,t) is the operator defined in Eq.(4.6). 

We now turn to S^°^^{ti,t2,t3]t) and remember [21] that the creation potential reads 
as = J^B'^i' ■ 'Ci^CL + ■ ■ - b^b); this shows that not only V^Jjji') = 0, which is the key 
property of the creation potential [21], but also {v\V^ = 0. The state {v\Ut-V^ is also 
equal to zero, but for a less trivial reason: This state is a zero-pair state since Ur destroys 
one pair while creates a pair. We can then insert \v){v\ in front of it to see that it 
reduces to zero since V^^l^) — 0- This shows that the non-zero terms of ta, ^3; 

must be constructed on {v\UrUr'V^Ul„\v) . By collecting them, we are left with three 
terms only that we split as 

t2, h; t) = t2, h; t) + ^2, h; t) , (4.9) 

where the first part is very similar to the two terms of S^^^\ti,t2,t3;t), 

S'„':°''\h,t2,h;t) = {v\Ut,e-'''-^''~'^^Ut,Vi^^^ + (^2^*3), (4.10) 
the second part reading 

Sy''\tut2,t3]t) = -(i;|f/t3e-^^-(*^-*2)f/t,e-^^-(*2-*)v;J,(t-ti)f/^|t;)e-*^™(*i-*) . (4.11) 

To get 7;;^^'^^^^) and 7;;^°"^(t), we are left with calculating S^^^'K^i, ^2, ^3; t) and Sg°''\ti, t2, ts; t) 
and integrating these quantities over (^3,^2,^1)- 

4.3 Calculation of 7^,^^^^(t) 

By using R^'^^^ given in Eq.(4.6), the expression of the semiconductor-photon coupling 
given in Eqs. (2.5,6) and Eq.(2.1), we can rewrite S^^^\ given in Eq.(4.8), as 

ji,j2,j3 ni,n2,n3 
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X {v\Bn,Dn,„,Bljv)F{ti-t,t2-ti,t-i-t2) , (4.12) 

where F{ti,T2,T3) reads as 

(4.13) 

with = ujj — En being the detuning of photon j with respect to exciton n. 

The matrix element in the RHS of Eq.(4.12) readily follows from one of the key equa- 
tions of the composite-boson many-body theory [20,21], namely, 

[^™. 5]] = E & i) + ("^ - ^)] Bi , (4.14) 

n 

where A is the Pauli scattering between two excitons. Since {v 

while Dmi\v) = 0, we get 

{v\Bn,Dn,rnBU\v) = A + {n, ^ n^) . (4.15) 

By inserting Eqs.(4.13,15) into Eq.(4.12), integrating over ti, T2, between — oo and 
and ultimately taking the complex conjugate of the result, we find the following compact 
expression for 7m^"^'(t) in terms of Pauli scatterings between two excitons, 

{3e-iujj ~iu)j +iuij )t , , (jl)* , , 02)* , , (is) 



X 



A (:;^ + (ni n^)] . (4.i6) 



4.4 Calculation of 7S°^^(i) 



(i) Let us start with given by Eq.(4.10). The procedure used for S*^ leads us to 

write 

rti rt2 , nUl) ,iiJ2) ii{j3)*p{3e+iujj-^+iujj^-iujj^)t 

/ dt2 rft3^:^-'(t„t„t3;t)= E E "fri^o v.-,o ^ 

J-00 J-00 Jl,j2,i3 '^l'"2,n3 + *^^i2"2Jl^ *^'i3»^3/ 

X G'(ti - t) , (4.17) 

with 

G'in) = e(3^+^'^n+^'^.2-*^.3n3)-i {v\Br,,Bn,Vi{Ti)BlJv) . (4.18) 

To go further, we use Eq.(3.6) for V^{ti) and, in its last factor, replace Hgc by En^. This 
allows us to formally perform the x integration as 

p~iHscT\ p—i{En^+Em)T\ 

Vi{n)Bl\v) = _ _ ViBl\v) . (4.19) 

^ sc ^n-j 
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Integration of G'{ti) over ri then yields 

dnG'in) = ] _ (.|i?AT-K^<|.) , (4.20) 

where the operator T' reduces to 

1 

X' = . (4-21) 

3e + iuj-^ + iuj^ — i^j.^n^ — iHgc 

(ii) We now turn to 5*^^°"' defined in Eq.(4.11). We again replace the t/'s by their 
expression (2.5,6), use Eq.(2.1) and integrate over ^3 and ^2- We find that this integral 
reads as Eq.(4.17), with G'{ti — t) replaced by G"{ti — t) given by 

Wn.B,.. , , . , , . JTT Vi{-T{)Bl\v) . (4.22) 



By inserting Eq.(4.19), with ri replaced by (-ti), into the above equation, we get an 
expression of G"{ti) which again can be formally integrated over ri. The result just reads 
like Eq.(4.20) with the operator T' replaced by T" given by 

T" = ^ ~ ^^^'^"3 _ 23) 

(3e + iojj^ + iojj^ — ifij^ns, — iHgc) (2e + iuj-^ + iuj^ — iHgc) 

(iii) The last step is to add the two contributions to 7i^°'^^(t) coming from the triple 
integrals over times of 3'^^°^^ and 5*^*"°*^^. We first note that 

T' + T" = : : . (4.24) 

2e + itUj-^ + iuj^ — iHgc 

By using Eq.(3.10) and taking the complex conjugate of the final result, we end with a 
nicely compact form of 7m°"^(^) iii view of the rather complicated calculation we had to 
perform to get it, namely 

J3e-iujj -iuij +iujj )t (ii)* (i2)* , , (is) 
^Coulj^^-j ^ ^ ^ Pni Pn2 f^ns 

ji,j2,j3 "I'^ains i^jim + ^j2m " ^jsm + 3ze)(r2j2n2 + ^f)(^i3"-3 ^ 

X y e f"' {v\Bn'B:^: i Bl Bl \v) . (4.25) 

5 Discussion 

When compared to the Pauli contribution given in Eq.(4.16), we see that this Coulomb 
contribution has a very similar structure. This leads us to write 7m(^) as 

- i (t-Jji -ojjg )i ,,{jl)* ,,{32)* ..(j'i) 
ry^^^{t) = e""^ y y - Pni Pn2 ^^n■i 

jij-ij-i "•i"2"3 (^iim ~^ ^j2m ~ ^jsm + S^c) (f2j2„,2 + ie) (fijg^g — ie) 

TZ^''''{nun2,ns)+TZ°''\ju32;nun2,ns)] , (5.1) 
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where F^^^'' and r^°"^ are physically linked to interactions with the exciton m, through 
Pauli exclusion and Coulomb interaction respectively. Their precise values are given by 

r^-V, n2, ns) = - [a -) + in, ^ n^)] , (5.2) 
5.1 Large detuning expansion of ^m{t) 

As previously said, the Pauli part of 7m(^) is fully compact since it only reads in terms of 
the Pauli scattering A In contrast, the Coulomb part cannot be obtained exactly 

because the full spectrum of the two-pair eigenstates is unknown. We can however note 
that this Coulomb term contains one additional ratio of Coulomb scattering divided by 
photon detuning. This makes it small far from resonance, i.e., when the concept of 
nonlinear susceptibility is relevant compared to absorption controlled by the Fermi golden 
rule, i.e., the poles of ■jmit)- 

In order to get the dominant contribution in the large detuning limit of r^°"'(ji, ni,n2, n^), 
we can make Hgc act on the right in Eq.(5.3) and, according to Eq.(3.4), replace it by 
En-^ + En2- By using Eq.(3.8), this leads to 

r^°"'(ji,j2;ni,n2,n3) ~ ^ ^- ^ \ , (5.4) 

where ^"^ is the "in" Coulomb exchange scattering of the composite-boson many- 

body theory [21]. More generally, Eqs.(3.4), (3.8) and (3.9) allow us to expand T^a^^Ui, 32] n,, n2, n^) 
in ratios of Coulomb scattering divided by photon detuning. 



5.2 Shiva diagram representation 

To better grasp the physics of this third order response to a multifrequency photon field, 
Shiva diagrams [21,22] turn out to be once more quite valuable. 

In a linear response, one photon with frequency Uj^ is coupled to the exciton m. 
This corresponds to the diagram of Fig. 1(a). Since with one photon, we can only have 
one exciton, the linear response shown in this diagram makes appear fjf:^^* / [uj — Em)'- 
The numerator /i^^-** comes from the transformation of one photon with frequency ujj^ 
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(we here drop its t dependence to simplify the notations) into the exciton with energy 
Em, while the denominator, necessary to have a dimensionless quantity, is, as usual, the 
difference between the energies before and after the transition, as shown by the vertical 
lines in Fig. 1(a); before the transition, we only have the photon uj-^, and after, we have 
the exciton m. 

If other photons are present, the excitons to which these photons are coupled can 
interact. The simplest interaction is a bare fermion exchange (see Fig. 1(b)). In this dia- 
gram, we have two couplings with the "in" photons cuj^ and toj^ which induce n^^"^* and 
l^n2^*- We also have one coupling with the "out" photon tOj^ which induces We 
obviously also have the Pauli scattering for fermion exchange A with (rii,n2) pos- 

sibly exchanged. The energy difference between the "in" state with photons cUj^ , ujj2 and 
the "out" state with exciton m and photon uj^ has to appear in the denominator. More- 
over, in order to have a dimensionless quantity, we need two other energy denominators 
since the /x's are energylike quantities. These denominators can only be the detunings 
associated to the excitations of the virtual excitons n2 and which can be seen as the 
energy difference between states before and after these couplings (small vertical lines in 
Fig. 1(b)). Consequently, the contribution of the diagram shown in Fig. 1(b) reads 

nih)* nih)* iiUa) r / X i 

the minus sign being standard when the number of exchanges is odd. This exactly is the 
result given in Eq.(4.16). 

In addition to this pure exchange, excitons can also interact through the Coulomb 
potential. If we restrict to one Coulomb process, we are led to diagrams (c) and (d) 
of Fig.l. These diagrams make appear the direct Coulomb scattering ^ and the 

exchange Coulomb scattering {^^ JJ^^ with a minus sign since in this last scattering, 
one exchange is present. However, as these scatterings are energylike quantities, we need 
one more energy denominator which turns out to be the energy difference between the 
"in" state with the two photons {uoj^^uoj^) and the "out" state with the two excitons ni 
and n2- Consequently, the contribution of diagrams (c) and (d) of Fig.l reads 



nih)* nih)* iiUs) 
h^ni 



^ji + ~ — En2 



(5.6) 
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We thus recover the expression of the large detuning hmit of 7m°"'(^) given by Eqs.(5.1) 
and (5.4). Other Coulomb interactions between the two excitons can be included along 
the same line to get the detuning expansion of the third order response function. 

We see that the above discussion based on Shiva diagrams for exciton interactions 
[21,22] makes the various terms of •jmit) rather obvious. 

5.3 Ground state exciton approximation 

From Eqs.(2.8) and (5.1), we see that (J^'^''(r))t contains a summation over excitons m, ni, 
77-2, ?T-3, with n = (z/„, Q,„) where z/„ is the relative motion index and Q„ the center-of-mass 
momentum of exciton n. 

(i) Since through the exciton-photon couplings relevant excitons {ni,n2,ns) have 
a momentum equal to the one of the photon to which they are coupled , while both Pauli 
and Coulomb scatterings conserve momentum, all exciton momenta in Eq.(5.1) are of the 
order of photon momenta, i.e., small on the characteristic scale of excitons. Consequently, 
we can take and all Qn's equal to zero. 

(ii) J(r) a priori contains contributions from all excitonic levels. For practical use, it 
is however possible to restrict the summation to the ground state exciton: The photons 
at hand usually have a detuning with respect to this ground state exciton which is rather 
large to avoid absorption on the sides of the exciton line, but not that large to possibly mix 
the various exciton levels in the nonlinear susceptibility. Moreover, due to the (r = 0|z/) 
factor appearing in j„^(r) and in these couplings are the largest when corresponds 
to the relative motion ground state {uq). For these reasons, it is possible to replace 
and all u^s by z/q, with one exception in the Pauli term as we now show. 

(iii) As a matter of fact, by looking at this Pauli term, Eq.(4.16), we see that the sum 
over z/„^ can be readily done, which is more valuable than replacing by uq. The only 
factors which depend on are the coupling n^j^^^* proportional to (i^nilr = 0) and the 
Pauli scattering A which, when all exciton momenta are equal to zero, reduces to 
[20,21] 

A (:;^ :i) = E(^™|k)(^^njk)(k|z/„,)(k|z.„,) , (5.7) 

k 

where (k|z/„) is the relative motion wave function of state n in momentum space. The 
closure relation J2,/rii ~ ^ then allows us to readily calculate the sum over z/„^: 
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After replacing n2, m and m by o = (i/q, Q = 0), we get Em /^^ ^*-^ (o nj = l^o^'^*-^o, 
where Ao is given by 

Ao = El(^o|k)P(z/o|k)j^f^. (5.8) 

On the other hand, in the large detuning limit of the Coulomb term given by Eqs.(5.1) 
and (5.4), we can replace m and the n's by o = {uq, Q = 0). Since (o o) = Pl]? we are 
left with the Coulomb exchange term in which appears ^™ (° °) given by [21] 

r (^ o) = 2E^k-k'|(^o|k)P [|(z/o|k')p - (z/o|k)(kVo)] . (5.9) 

k,k' 

(iv) The third order current density can thus be approximated by 



(G.A^J* (G.A^J* (G.A,3, 



X V"— jiy V"— j2y V"-— j3y _^ ^^^^ (5.10) 

(r^j^o + - 1)^30 + 3ze)(^20 + «e)(^30 - ze) 

where the prefactor Kq reads as 

Ko = 2G (a, + t^^tFWI = . (5.11) 

By replacing the ground state wave function in momentum space (k|z/o) and the Fourier 
transform of the Coulomb potential Vq by their values in 3D or in 2D, in Eqs.(5.8) and 
(5.9), we can easily show that Aq = Xoicix/ L)^ and ^™ (° °) = —^dRx{clx/ L)^ , where 
and Ox are the 3D exciton Rydberg and Bohr radius, the prefactors Xd and being equal 
to 77r/2 and 267r/3 in 3D, while they are equal to 27r/7 and (87r - 3157rV512) ~ 6.1 in 2D. 
Note that, due to the factor appearing in these scatterings, Kq is indeed volume 

free. 



6 Conclusion 

This paper shows rather in details how to calculate the current density induced by a 
multifrequency photon field, starting from its standard textbook expression as a third 
order time integral of a triple commutator. It ends by a rather simple approximate 
expression of this current density when the ground state exciton plays a dominant role. 

The procedure which is here reported heavily relies on the tools we have developed 
in the many-body theory for composite excitons [20,21]. These tools mainly allow us to 

18 



extract the volume linear terms appearing in the standard expression of this nonlinear 
susceptibility and to show their exact cancellation. This exact cancellation, which is 
a necessary requirement since susceptibilities are intensive quantities, stayed an open 
problem for decades. The current density ends by reading in terms of Pauli scatterings for 
carrier exchange between two excitons in the absence of carrier interaction, and expands in 
terms of direct and exchange Coulomb scatterings describing Coulomb interaction between 
the carriers of two excitons. Once again, Shiva diagrams greatly help to understand the 
various contributions to the third order response. 

Although the calculation of this third order susceptibility is conceptually rather trivial 
because two excitons only are involved, so that it is not at all a A^-body problem, the 
precise algebra necessary to combine all the terms to get the volume-free final result in 
a compact form is not that straightforward. This is why we have found useful to report 
this calculation rather in details for the reader to possibly reproduce it without difficulty. 
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APPENDICES 



A Current density and semiconductor-photon inter- 
action 

Although most authors describe nonhnear optical effects in terms of polarization, it is 
of importance to note that the physically relevant quantity for solid state physics is not 
the polarization, but the current density. Indeed, as explained in details for example 
in Cohen-Tannoudji et al.'s textbooks [32], the change from A.p to E.d relies on an 
integration by parts which is only valid for wave functions cancelling at infinity, as in the 
case of atomic systems, but not for solids with lattice periodicity. Since we could not find 
a textbook in which the following derivation is performed from scratch in a simple way, 
it appeared to us as useful to have it printed somewhere since the photoinduced current 
density is the key for any clean description of optical nonlinear effects in semiconductors. 

A.l Current density 

The current density of a set of classical particles located at r„ with charge q = — |e| and 
velocities v„, is given by 

J(r) =5]gv„5(r-r„) . (A.l) 

n 

Its Fourier transform thus reads 

Jq= /t/rJ(r)e^Q-"- = 5]gv„e^Q-'-" . (A.2) 

For quantum particles, the associated operator in first quantization is obtained, as usual, 
by symmetrizing the classical expression with respect to (r,p); since v„ = p„,/m, this 
operator thus reads 

Jq = ^ E (Pn e^"^-^" + e^Q-- p„) = JIq . (A.3) 

To write it in second quantization, we use the valence and conduction Bloch states creation 
operators d\ where n = (c, k) or {v, k), with k quantized in 27r/L for a periodic solid of 
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size L. The standard second quantization procedure then leads to 

Jq = J2iQi^'^n)(^l'dn , (A.4) 

n',n 

where the prefactors are given by 

jQ(n',n) = J dv<Pl,{v)^ (pe^^^ + e^^^p) 0,(r) . (A.5) 

Interband processes correspond either to n' = (f,k') and n = (c, k), or n' = (c, k') 
and n = (f , k). The wave function 0„(r) for Bloch state n = (c, k), with r = R + u and u 
restricted to a unit cell, R running over all unit cells, reads 0ck(i") = iv~^/^e''''*-^''~"Vck(u). 
By using a similar form for n' = [v, k'), integration over R in Eq.(A.5) leads to k' = k+Q. 

If we only keep these interband terms, which are the ones of interest for optical phe- 
nomena as they can lead to resonant processes with photons, we get 

Jq = E ^(k)4,k+Q4,k + E j?^(k)<k4,k-Q . (A.6) 

k k 

Since Jq = jLq due to Eq.(A.3), the prefactors in J'^^^ must be such that jQ"'*(k) = 
j^Q(k) . We then note that (i) the Q's of interest are usually small, (ii) the operator 
p has nonzero matrix elements between valence and conduction states, since these states 
have different parities, the operator p being odd. This allows us to neglect the (Q, k) 
dependence of jQ'^^(k) and to replace it by its finite value limit obtained for Q = = k. 
Equation (A.6) can then be replaced by 

Jitter ^ _G ^ dU+Q4,k -G*Y: 4,k4,k-Q , (A.7) 
k k 

where G, often called Kane vector [30], is given by 

G= / duv7,o(u)(-zV)v9eo(u) , (A.8) 

with Uc being a unit cell. 

The next step is to turn from valence-conduction electrons to electrons and holes 
according to dl ^^ = aj^ and c?]^ ^ = b-^. It is then appropriate to note that 

E 4,kC^f,k-Q = E «k^-k+Q = E "p+aeQ^-p+ahQ ' (A-9) 
k k P 
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with = 1 — a/i = me/{me + rrih). By using the hnk between free pairs and excitons, 
namely, 

V 

and by noting that (r|p) = e*P'''L~'^/^, so that 

E(Hp) = E(HP)(p|r = 0)L^/2 = (Hr = 0)^^/^ , (A.U) 
p p 

the sum over k in Eq.(A.9) ends by reading in terms of exciton operators as 

E<kC?.,k-Q = E^™^Q™,Q(^-Ir = 0)i:^/' . (A.12) 

k m 

Consequently, the interband part of the current density Fourier transform Jq can be 
rewritten as 

jinter ^ _ Q J] B„5q^,_q (r = ^\v„^L^I^ - G* E ^^^Q.,Q (^^1^ = 0)L^/^ . (A.13) 

m m 

Note that these two terms are not complex conjugate, in agreement with Jq = JIq, the 
Fourier transform of the current density being not real, as seen from Eqs.(A.2) and (A. 3). 
If we now come back to the current density in r space, its operator, given by 

has an interband contribution which ultimately reads in terms of exciton operators as 

j'°*''(r) = Ej-(r)i?m + h.c. , (A.15) 

m 

in which the prefactors are given by 

j^(r) = -Ge'Q-'-(r = 0|i.^)L"^/2 . (A.16) 

We can check that the current density has a classical expression, Eq.(A.l), which is real, 
so that it is represented by a Hermitian operator, in contrast with its Fourier transform 
Jq. 

A. 2 Link with semiconductor-photon interaction 

The kinetic part of a set of electrons in an electromagnetic field with vector potential 
A(r, t) is known to be given by 

^0 = E (Pn - 9A(r„, t)f /2m , (A.17) 

n 
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so that the electron-photon interaction in first quantization appears as 

= E (Pn-A(r„, t) + A(r„, t).p„) . (A.18) 

Since the replacement of A.p by E.d [32] cannot be done for the extended wave func- 
tions of solid state physics, because it relies on integration by parts with exact cancellation 
of the integrated term, the only correct way to describe semiconductor-photon interaction 
is to stay with A.p, as we now do. However, it can be of interest to note that the symme- 
try properties of p and r are in fact the same. It turns out that these symmetry properties 
are the only ones really used for most optical effects, the size of the coupling being more 
or less an adjustable parameter. This is why results obtained from E.d instead of A.p 
can look as correct, although derived from an incorrect first line. Nevertheless, since there 
is no specific advantage of working with E.d instead of A.p, we do not see the interest 
of using E.d in semiconductor physics, except for some readers who may find the dipole 
density physically more "appealing" than the current density. 

The vector potential A(r,t) for a set of photon fields with amplitude |Aj| made of 
{u!j,Clj) photons, introduced adiabatically from t = —oo on the same scale 1/e for sim- 
plicity, reads as 

A(r, t) = e'' (A,e^("^*-Q^-'-) + c.c.) , (A.19) 
j 

so that the electron-photon interaction splits as Wt = J2j W}''^ , with 

= _eet^icu,tj^^_J_ rp^e-iQ..r„ ^ e-^Q.-rnp^A ^ Yi.c. . (A.20) 
2m „ 

When compared to Eq.(A.3), the above equation readily shows that the electron-photon 
interaction is simply related to the Fourier transform of the current density through 

Equation (A. 12) for the interband contribution to Jq, then gives this coupling in second 
quantization as 

^ ^et^iu^.t (^jj, ^ rj.]\^ ^ Yi.c. , (A.22) 

where we have set 

U, = E f4^Brn T^ = Y: r^Bl , (A.23) 

m m 
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the prefactors being related to the Kane vector G through 

f4l> = G.A, <5q,„,q^ (r = 01^.^) , (A.24) 

= G*.A, 5q„_q^. {um\r = 0) . (A.25) 

Among these two sets of terms, only {Uj, t/j) can lead to possible resonance with photons. 
This is why, for photons close to the band gap, we are led to perform the so-called "rotating 
wave approximation" [31] which reduces the semiconductor-photon interaction to 

^0-) ^ e^'e'^^^Vj + h.c. . (A.26) 
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B Commutator expansion of operator mean value 

Let us consider an operator At, which possibly depends on time, and a state \ipt) which 
obeys the Shrodinger equation, id\ipt)/dt = {Hq + Wt)\4't)- We want to calculate the 
mean value {At)t = {ipt\At\ipt) as an expansion in powers of Wt- 

A convenient way to do it is to note that this mean value, which also reads {ipt\e~^^°^ Ate^^^^lipt) , 
where Z = ^^Hot^^-iHot ^y^q Heisenberg representation of operator Z, can be written as 

m m 

SO that it also reads 

{At)t = TT{ARt) , (B.2) 

where Rt = \ipt){'ipt\ is the projector over state \ipt)- 

The expansion of Rt in powers of Wt is easy to derive by noting that d/dtRt = 
—i[Wt,Rt], so that if we expand Rt in powers of Wt as J2n=oRt^\ the Rt^^'s are related 



by = -z[Wt,R 



(n-l)n 



For Wt: 



and Rt= 



R, this set of differential 



equations leads to 



R 



(0) 



R 



R 



(1) 



dh[Wt„R] 



R 



(2) 



dti 



dU 



Wt,,[Wt„R] 



(B.3) 



and so on. . 



This allows us to expand {At) t as E+ro(A)^\ where {At)r' = Tr(Ai?rO- This trace 
is easy to write in a compact form by noting that TtAB = AijBji = TiBA, so that 
Tt{AB)C = TtC{AB), which leads to Ti A[B, C] = Tt[A, B]C, and so on. Consequently, 

{At)i'' 



An) 



dhTT[At[Wt,,R_ 
dt.Ti ([At,Wt,]R 



(B.4) 

If = \v), with |t;) being Hq eigenstate, we do have \v){v\ = e^^°^\v){v\e~^^'^^ = R, 

so that Tr ([At, Wt]Rj = {v\[At, Wt]\v). A similar procedure for higher order terms leads 

to write the third order term in Wt as 

rt rti (-12 



{At) 



(3) 



-I 



dti 



dto 



dt3{v\ At,Wt, 



Wi 



t2 



Wi 



\v). 



(B.5) 
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Figure 1: (a) Linear response to a photon field: One photon with frequency Uj-^ transforms 
into one exciton in state m. (b) Third order response: Two photons uj^ and Uj^ transform 
into a third photon uj.^ and one exciton in state m. This can be done by just carrier 
exchange between the virtual excitons to which photons {ujj-^,ujj^) are coupled. (c,d) This 
third order response can also contain Coulomb processes. The one with one Coulomb 
interaction and no exchange is shown in diagram (c), while the one with one Coulomb 
interaction and a carrier exchange corresponds to diagram (d). 
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